Photoassociation of sodium in a Bose-Einstein condensate 
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Abstract 

We report on the formation of ultra-cold Na2 molecules using single-photon photoassociation 
of a Bose-Einstein condensate. The photoassociation rate, linewidth and light shift of the J = 1, 
V = 135 vibrational level of the A^S^ molecular bound state have been measured. We find that 
the photoassociation rate constant increases linearly with intensity, even where it is predicted that 
many-body effects might limit the rate. Our observations are everywhere in good agreement with 
a two-body theory having no free parameters. 

PACS numbers: 03.75.Fi, 33.70.Jg, 34.20.Cf, 33.20.Kf 
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Bose-Einstein condensates (BECs) of atomic gases are versatile systems for the 
experimental study of quantum behavior. Of particular interest are the suggestions for the 
coherent coupling of a BEG of atoms with a BEG of molecules ||, |[ and the possibility 
of creating entangled pairs of atoms in the BEG via coupling with molecular levels [0]. Two- 
photon photoassociation processes using stimulated Raman transitions have formed ground 
state molecules from ground state atoms in a BEG Q but at very low rates. Here we 
explore the fundamental upper limits of molecule formation by making them at high rates 
using the elementary process of single-photon photoassociation. 

In single-photon photoassociation two atoms collide in the presence of a light field and 
form an excited state molecule. Photoassociative spectroscopy is used extensively to study 



collisions between laser cooled atoms [T^. Photoassociation in a BEG presents quite a 
different regime: the collision energies are orders of magnitude lower than in a laser cooled 
sample (the de Broglie wavelength is as big as the sample) and the densities are higher. This 
puts us in a regime where many-body effects may be important. 

We concentrate on a particular photoassociation transition and measure the photoas- 
sociation spectra for various intensities and durations of the light pulse. From these we 
determine the photoassociation rate, lineshape and the shift of the resonance. Finally we 
examine various limits on the photoassociation rate. 

Figure ^ shows the photoassociation process. The molecular level chosen for study is the 
J = 1, V = 135, rotational- vibrational level of the A^S^ Na2 molecular state, excited from 



free atoms by a laser frequency of 16913.37(2) cm^^ We chose this level because its 

detuning from the Di resonance is far enough (43 cm~^) for atomic absorption to be negligible 
and because earlier experiments using laser cooled atoms indicated a high photoassociation 
rate. The lifetime of our excited molecules is about 8.6 ns. Decay of excited molecules into 
hot atoms or ground state molecules constitutes loss from the condensate. This loss is how 
we detect photoassociation. 

We prepare an almost pure condensate of ~ 4 x 10^ sodium atoms in the \F = l,mp = 
— 1) ground state with a peak density of rio ~ 4 x 10^^ cm"'^. The condensate is held in an 
anisotropic magnetic TOP [|11|] trap with oscillation frequencies of lUx/\/2 = ujy = \/2ujz = 
27r X 198 Hz and corresponding Thomas- Fermi radii of \/2rx = ry = /V2 = 15 fim. 

To induce photoassociation we illuminate the condensate with a Gaussian laser beam 
focused to 120 fim FWHM at the condensate. The peak intensity is varied from 50 to 1200 
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W-cm~^. The polarization is linear and parallel with the rotation (z) axis of the TOP trap 
bias field. The hght is applied as a square pulse for between 10 and 400 /xs with rise and fall 
times of less than 0.5 /is. 

The condensate number is measured using phase contrast imaging [|T2] , taking two images 
before and two images after the photoassociation pulse to determine loss. The imaging occurs 
at 40 ms intervals using a 100 /is probe pulse from a laser tuned 1.78 GHz to the red of 
the 3Si/2,-F = 1 ^ 3P3/2,-F = 2 transition. The imaging rate is limited by the readout 
time of the CCD camera. The photoassociation pulse begins halfway between the second 
and third images. We use multiple imaging pulses to improve statistics and to partially 
correct for small losses other than those due to the photoassociation pulse. These losses are 



typically 4% between images. Once the number of atoms is extracted from the images ||23 
we calculate /, the fraction of atoms remaining after the photoassociation pulse. 

Figure shows a typical photoassociation spectrum. Each point represents a freshly 
prepared condensate. We use a Fabry-Perot etalon and a reference laser locked to an atomic 
Na line to measure differences in the photoassociation frequency with a precision of 5 MHz. 
The laser linewidth is < 3 MHz. All detunings quoted are relative to the center of the 
photoassociation line in the low intensity limit. For small trap loss we expect the line to be 
Lorentzian (in contrast to photoassociation lines in an uncondensed thermal sample where 
the kinetic energy distribution distorts the line shape [T^). For significant trap loss, as in 
figure 1^, one must account for the change of the density profile during the photoassociation 
pulse. 

A two-body collisional loss process changes the local atomic density as: h{t, r) = 
—K{I,u)n^{t,r), where K{I,lj) is the intensity and frequency dependent photoassocia- 
tion rate constant. Because the characteristic time for the motion of the atoms, the trap 
oscillation period, is long compared to the photoassociation pulse, we assume that the local 
density only changes due to photoassociation. We then obtain: 

"(''^) = l + ir(/,..MO,r)f 

The density distribution flattens with time. Spatially integrating equation (|I|), assuming an 
initial, parabolic (Thomas- Fermi), density distribution and a uniform intensity J, leads to 
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an expression for the fraction of atoms remaining in the condensate 

-(l + r/)i/^anh-i(vW+^)}, (2) 

where r/ = K{I ,uj)n{0,O)t = K„,{I)n{0 , 0)t / {1 + {2{uj - uJo{I))/-i{I)f). We use a three 
parameter fit of equation (Q) to the spectra to extract the on-resonance rate constant Km{I), 
effective hne width 7(1) and central frequency ti;o(/) (for example, the dotted line in figure 0). 
The fit is good. To further verify equation (0) we plot the measured 1 — / as a function of 
pulse length for / = 140 W-cm~^ and u = Uq, along with a one parameter (Km) fit to the 
data (figure The error bars are the fitting uncertainty. 

By fitting spectra obtained at various intensities we measure Km{I), and ti;o(/). We 
calculate the unbroadened molecular linewidth of the chosen state to be 7o/27r = 18.5 MHz 
(nearly twice the atomic linewidth). This is in good agreement with the measured linewidth 
of 19.5(25) MHz in the low intensity limit where it is independent of intensity. At higher 
intensities we observe broadening with a maximum linewidth of 60 MHz at 1 kW-cm~^. 
Homogeneous power broadening is calculated to be three orders of magnitude too low to 
explain this width. It is, however, partially explained by differential light shifts of the 
unresolved molecular hyperfine states. These states are calculated to be split by less than 
1 MHz at low intensities and about 30 MHz at our maximum intensity. Another possible 
contribution is the inhomogeneity of the photoassociation beam intensity combined with 
the large light shift (discussed below). Variations due to either local spatial inhomogeneity 
(e.g. interference fringes) or displacement of the cloud from the center of the Gaussian beam 
could account for the extra width. Assuming these inhomogeneous broadening mechanisms 
do not change the area of the line (verified by a simulation), we take the on-resonance 
photoassociation rate constant to be Kq{I) = Km{I)l{I) 0). 

Figure § shows corrected and uncorrected as a function of intensity (for vari- 
ous pulse lengths). The error bars are the fitting uncertainties. Once we correct for 
the inhomogeneous broadening we get a linear dependence on intensity with a slope of 
dKo/dl = 3.5(2)(10) x 10~^° (cm^ ■ s~^)/(kW ■ cm~^). The uncertainties are due, respec- 
tively, to fitting and to the combined uncertainties in the measurement of the intensity 
and density. For intensities above 1.2 kW-cm~^, which we could only achieve by more 
tightly focusing the photoassociation laser, atomic dipole forces significantly perturb the 
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condensate, thwarting meaningful measurements. A coupled-channels, two-body scattering 
calculation with no adjustable parameters yields a photoassociation rate constant of 
dK^/dl = 4.1 X 10^^'^ (cm'^ ■ s^^)/(kW ■ cm^^) for our range of intensities. This includes a 
factor of 2 decrease relative to K for a non-condensed gas and agrees well with our experi- 
mental result. 

We study the photoassociation light shift, previously observed in a non-condensed gas 
jTSl , in a set of experiments where the total fiuence (intensity x pulse length) of the pulse 
was kept constant, to maintain the depth of the photoassociation dip in an easily observable 
regime. The results are shown in figure ^. The measured light shift is —164(35) MHz /^^, 
which leads at high intensity to a shift large compared to the linewidth. The principal 
contribution to the uncertainty is the intensity calibration. During preparation of this work 



we became aware of similar results in lithium 16 



While the strength of the photoassociation resonance is dominated by s-wave scattering, 
the dominant contribution to the light shift is due to a d-wave shape resonance. A the- 
oretical calculation of the light shifts using equation (3.7) of reference ^7^, including the 
effect of the d-wave shape resonance embedded in the continuum, gives an average value of 
— 130 MHz with a spread of ±13 MHz due to the hyperfine structure. 

We now consider the upper limit to the photoassociation rate constant K [this implies a 
lower limit on the photoassociation time r = {Kn)^^. If one uses a semi-classical theory 
that is commonly applied to collisions of laser-cooled atoms |TB| then Kq = crv = tcR^Pv, 



where Rc is the Condon radius (see figure |lD and P is the probability of a photoassociative 
transition with a maximum value of 1. If we take the relative velocity v of the atoms to be 
h/(2mry) = 0.6 mm-s^^, where m is the atomic mass, then the maximum photoassociation 
rate constant is four orders of magnitude lower than our highest measured value. This reveals 
the inadequacy of a semi-classical approach, which fails to take into account threshold laws 

Quantum theories for the photoassociation rate constant can be compared by expressing 
K as K = {h/m)L, where L is a characteristic length. Two-body s-wave scattering theory 
for a BEC gives Lg = \S{k)\'^/k where hk is the relative collision momentum and S{k) is 
the S'-matrix element for atom loss via photoassociation. References [O, Il9[ show that, on 
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resonance, 



where hr{k,I) = 2TT\{e\hQ\k)\'^ is the Fermi-golden-rule stimulated-decay width of the ex- 
cited molecular state |e) due to the optical coupling hil oc vT with the colliding atoms. 
Since F oc A; as A; ^ 0, and F/70 < 0.001 in our range of power and collision energy, Lg is 
independent of k. We calculate 7o/27r = 18.5 MHz. Lg is linear in / for our experimental 
conditions and dLs/dl is calculated to be 24 nm/ This gives the above-quoted rate con- 
stant in good agreement with the experiment. In our power range, can be significantly 
larger than the Condon point for the transition, 2.0 nm. Note that equation (|^) shows that 
Ls{I) will saturate with increasing I and decrease for sufficiently large /. 

The upper limit to the two-body quantum K is the unitarity limit where \S\'^ = 1 so 
Lu = 1/k = X/27T, where A is the de Broglie wavelength. Since A is on the order of the BEC 
size Ls/Lu ^ 1 so that our experiment is well below the unitarity constraint. 

Recent many-body theoretical work 21] has suggested an upper limit to in a 

BEC of ~ m-^J^ where Lj = ^2]!^ n~^^^ is the mean distance between particles. 
One might question if two-body scattering methods are applicable at densities where Lg 
becomes larger than Lj. At our maximum density Lj = 22 nm, so Lg/Lj ^ 1 at our 
highest intensities. This is the regime where one might expect two-body theory to start to 
fail. Nevertheless the linearity of K{I) at our highest intensities (figure ^ suggests that 
two-body theory continues to be valid. 

Larger values of Lg/Lj might be accessible by a modification to our experimental design. 
We can use the atomic dipole force (which currently limits our ability to use high intensities) 
to our advantage by trapping the atoms with the photoassociation laser. Without changing 
the atomic dipole forces, the laser can be suddenly brought from far off molecular resonance 
to on molecular resonance to induce photoassociation. Difficulties due to the molecular light 
shift might be reduced by finding a transition with a smaller light shift. 

In conclusion, we have measured the single-photon photoassociation in a BEC, in good 
agreement with two-body theory. This agreement represents a confirmation of the factor- 
of-two reduction for a two-body inelastic process in a BEC. The characteristic time for 
photoassociation is as short as 5 /is, much shorter than the 100 /is to traverse the mean 
distance between atoms, another demonstration of the extreme quantum nature of the col- 
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lisions. Our largest rate is still much smaller than the unitarity limit, but is about equal to 
a limit suggested on the basis of many-body effects; we do not see evidence that the rate 
saturates at this limit. 
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FIG. 1: The two-atom potentials for the ground state and the excited state used for photoassoci- 
ation. The atoms are initially unbound and on the ground state asymptote and are excited into 
the J = 1, V = 135 level. Prom there they decay and are lost from the condensate. Rc = 2.0 nm 
is the Condon radius, the internuclear separation where the energy of a resonant photon matches 
the difference between the potentials. 
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FIG. 2: A typical photoassociation loss spectrum. A 140 W-cm ^ pulse was applied for 100 fis. 
The dotted line is a fit to equation (0). The uncertainty in the frequency for each point is 5MHz. 
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FIG. 3: The dependence of the maximum condensate loss on the photoassociation pulse length for 
/ = 140 W-cm~^; the curve is a fit of equation (|2|). 
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FIG. 4: Photoassociation rate constant as a function of intensity. The corrected data has been 
adjusted to account for the inhomogeneous broadening. 
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FIG. 5: The light shift of the resonance as a function of laser intensity. 
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